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The Milky Way's dark matter halo is thought to contain large numbers of smaller subhalos. These 
objects can contain very high densities of dark matter, and produce potentially observable fluxes 
of gamma rays. In this article, we study the gamma ray sources in the Fermi Gamma Ray Space 
Telescope's recently published First Source Catalog, and attempt to determine whether this catalog 
might contain a population of dark matter subhalos. We find that, while approximately 20-60 of the 
catalog's unidentified sources could plausibly be dark matter subhalos, such a population cannot be 
clearly identified as such at this time. From the properties of the sources in the First Source Catalog, 
we derive limits on the dark matter's annihilation cross section that are comparably stringent to 
those derived from recent observations of dwarf spheroidal galaxies. 
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I. INTRODUCTION 

In the standard theory of hierarchical structure for- 
mation, cold dark matter particles gather to form small 
halos, which later merge to form ever more massive ha- 
los. As a result of this process, dark matter halos (such 
that that hosting the Milky Way) are expected to contain 
many subhalos, ranging from dwarf spheroidal galaxies 
(M > IO^Mq) down to objects with masses as small as 
10~^ — 10~^ Mq. The minimum mass of dark matter sub- 
halos is determined by the temperature at which the dark 
matter particles became kinetically decoupled from the 
cosmic neutrino background, which in turn depends on 
the particle physics of the dark matter candidate. 

With the exception of dwarf spheroidal galaxies, efforts 
to directly observe dark matter subhalos have failed to 
detect any such objects. A promising method to search 
for nearby dark matter substructures is to use gamma 
ray telescopes, such as the Fermi Gamma Ray Space 
Telescope (FGST), to detect the products of dark mat- 
ter annihilations. To such a telescope, relatively large 
(~ 10^ - 10"^ Mo) and nearby (- 0.01 - 10 kpc) subha- 
los could appear as as bright and point-like gamma ray 
sources, without counterparts in other wavelengths. For 
reasonable assumptions (based on the results of numeri- 
cal simulations) regarding the number of subhalos in the 
Milky Way and the dark matter distribution within those 
subhalos, one is led to expect a handful of subhalos to 
be observable by the FGST. For example, a dark matter 
candidate in the form of a 50 GeV particle with a annihi- 
lation cross section oi av ^ 3 x 10~^^ cm'^/s is predicted 
to provide on the order of a few subhalos that would be 
observable at or above the 5cr level [1]. 

The FGST collaboration has recently published a cat- 
alog of 1451 point sources 2 , including 630 that are not 
associated with sources in other astronomical catalogs. 
Of these unidentified sources, 368 have been detected 
with greater than 5a significance, and are more than 10 
degrees away from the Galactic Plane. In this article, we 
study this collection of gamma ray sources and consider 



whether any significant fraction might be the result of 
dark matter annihilations taking place in nearby subha- 
los. 

The remainder of this article is structured as follows. 
In Sec. HH we calculate the number of dark matter sub- 
halos predicted to be observable by FGST, for various 
astrophysical assumptions and particle dark matter can- 
didates. In Sec. Ill we determine which of the objects 



contained within the FGST's First Source Catalog can 
be fit by dark matter annihilations, for various combi- 
nations of mass and annihilation channel. We use this 



information in Sec. IV to derive upper limits on the dark 
matter's annihilation cross section. In Sec. |V) we discuss 
some of the more interesting features of the First Source 
Catalog, and consider whether they might result from a 
population of relatively large and nearby subhalos. We 
summarize our results in Sec. IVII 



II. 



GAMMA RAYS FROM NEARBY DARK 
MATTER SUBHALOS 



Following the results of numerical simulations, we 
begin by assuming that the Galactic Halo contains 
dark matter subhalos with a mass function given by 
dNn/dMh oc M^^, down to some minimum mass, nor- 
malized such that 10% of the total mass of the halo is 
found in IO^-IO^^Mq subhalos For a minimum sub- 
halo mass of one Earth mass, this normalization cor- 
responds to a total of ^ 5 X 10^^ subhalos within the 
Milky Way, which collectively make up about half of 
our Galaxy's total mass. In the local neighborhood, this 
implies a number density of approximately 34 (roughly 
Earth mass) subhalos per cubic parsec. The precise value 
for the mimimum subhalo mass depends on the character- 
istics of the dark matter particle [4]. As this study relies 
primarily on the properties of the largest mass subhalos, 
the minimum mass adopted has no significant impact on 
our results. 

We further assume that the subhalos are described by 
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the following density profile, as supported by the Via 
Lactea simulation ^5j: 
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where Rg is the scale radius of the subhalo. One should 
note, however, that the halo profiles describing the in- 
ner volumes of dark matter subhalos are difficult to re- 
solve using present simulations, and groups other than 
Via Lactea collaboration (notably, the Aquarius Project) 
have reported somewhat less steep inner profiles [6]. 
With this in mind, we will also calculate results using 
an Einasto et al. profile [7]: 



p(r) c>c exp 



\-- 




a 


kRs) _ 



(2) 



where a ~ 0.17. 

The concentration of a halo is defined as the ratio of 
its virial radius to its scale radius. To estimate the con- 
centration of a halo of a given mass, we use the ana- 
lytic model of Bullock et al. [5]. This model estimates 
a concentration of 21 (27) for a subhalo of mass 1O^M0 
(IO^Mq). For large halo masses, the results of numeri- 
cal simulations are in good agreement with this model. 
Smaller subhalos, currently beyond the resolution of such 
simulations, are largely irrelevant to our study. There are 
also indications from numerical simulations that subhalo 
in the inner volumes of their host halo tend to have, on 
average, higher concentrations [9 . To be conservative we 
do not include this effect in our calculations. 

One should keep in mind that considerable halo-to- 
halo variation in the concentration and shape of subhalo 
profiles has been observed in numerical simulations. The 
model of Bullock et al., for example, only provides a mea- 
sure of the average concentration of a subhalo of a given 
mass. The probability of a halo having a given concentra- 
tion can be modelled by a lognormal distribution, with 
a dispersion of Cc ~ 0.24 8 . Although we adopt mean 
concentration values in our calculations, we estimate that 
halo-to-halo fluctuations will increase the number of ob- 
servable subhalos by a factor of about two in most cases 
(see also Ref. [T]). Note that, throughout this study, ref- 
erences to subhalo masses denote the mass prior to loss 
through tidal stripping. 

Subhalos in the local volume of the Milky Way are 
likely to have had a large fraction of their mass stripped 
through tidal interactions with other halos and stars. 
This primarily impacts a subhalo's outer density pro- 
file, leaving its denser, more tightly bound inner cusp 
intact [in]. As the default assumption throughout our 
calculations, we assume that nearby subhalos have lost 
the outermost 99% of their total mass (although this only 
modestly impacts the overall annihilation rate). 

Just as the halo of the Milky Way contains many sub- 
halos, subhalos themselves are also expected to contain 
smaller bound dark matter structures within their vol- 
umes. This can lead to an overall "boost factor" to the 



dark matter annihilation rate from such objects. To be 
conservative, we do not include any such boost in our 
(default) calculations. 

The rate of gamma rays produced from dark matter 
annihilations taking place in a nearby subhalo is given 
by 



2'^DM 
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where av and moM are the annihilation cross section 
and mass of the dark matter particle, respectively, and 
the second integral is performed over the volume of the 
subhalo. dN^/dEry is the spectrum of gamma rays pro- 
duced per dark matter annihilation, which depends on 
the dominant annihilation channel(s) and on the mass 
of the dark matter particle (we use PYTHIA [T2^ as im- 
plemented in DARKSUSY [T3] to calculate the gamma 
ray spectrum). In the left frame of Fig. [T] we show the 
gamma ray luminosity from a subhalo, as a function of 
the subhalo's mass. Results are shown for 99%, 90% and 
0% mass loss. 

In order for the gamma ray annihilation products from 
a subhalo to constitute a source that could potentially 
appear in the FGST point source catalog, the subhalo 
must be both sufficiently bright and sufficiently compact 
to mimic a point source. To estimate the number of 
events from a subhalo observed by FGST, we multiply 
the gamma ray flux by an effective area of 6800 cm^, 
and a coverage of 20% of the sky at any given time. Al- 
though the detectability of a given gamma ray source 
depends somewhat on its spectral shape and its location 
in the sky, we can roughly estimate how bright a given 
subhalo must be to be detectable at high significance by 
FGST. In particular, the diffuse gamma ray flux mea- 
sured by FGST generates approximately 20 events per 
year per square degree above 1 GeV over galactic lati- 
tudes of 10° < |6| < 20°, and about 60-70 events per 
year per square degree above 1 GeV over galactic lati- 
tudes of |5| > 60°. In these two regions of sky, we esti- 
mate that 5-\/20 « 20 or 5-\/60 « 40 signal events per year 
above 1 GeV would be required in order for a subhalo to 
be potentially discovered with 5f7 significance. Based on 
this estimate, we conservatively classify any subhalo that 
produces more than 50 events above 1 GeV per year at 
FGST to be potentially detectable. In the right frame of 
Fig. [T] subhalos below the solid line are sufficiently bright 
to be potentially detected by FGST according to this cri- 
teria. Although the smallest subhalos (sometimes called 
"microhalos" ) are much more numerous, individual sub- 
solar mass halos are unlikely to provide observable fiuxes 
gamma rays (also see, for example, Ref. 15 ). We find 
that most of the subhalos potentially detectable by FGST 
have relatively large masses, M ~ 10^ — lO^M©. 

To appear as a point source to the FGST, the angular 
extent of the halo must not be much larger than the tele- 
scope's point spread function. To be considered point- 
like, we require that 95% of the photons from a sub- 
halo come from within a 2° radius (approximately the 
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FIG. 1: Left frame: The number of gamma rays per second produced from dark matter annihilations in a subhalo as a function 
of the subhalo's mass (prior to any mass loss). The solid and dashed lines correspond to the cases in which the outermost 99% 
or 90% of the original mass of the subhalo has been lost through tidal interactions, respectively. The dotted line (which may be 
difficult to distinguish from the dashed line) represents the zero mass loss case. Right frame: The range of halo masses (again, 
prior to any mass loss) and distances for which the gamma ray annihilation products from a subhalo may be detectable by 
FGST and appear approximately as a point source. To be considered potentially observable, we require more than 50 events 
per year above 1 GeV from a given subhalo. We have assumed that 99% of the original subhalo's mass has been lost. To be 
considered point-like, we require that more than 95% of the photons from a subhalo be concentrated within 2° (approximately 
the 95% containment angle for 1 GeV photons at FGST). In each frame, we have adopted a dark matter profile as described 
in Eq. [l] and considered a dark matter particle with ttidm = 50 GeV, cru = 3 x 10~^^ cm'^/s (the typical value for a thermal 
relic), and that annihilates primarily to bb. 



95% containment angle for 1 GeV photons at FGST). In 
Fig. [l] only those halos below the solid line and above 
the dashed line are both sufficiently bright and point-like 
to potentially appear in the FGST point source cata- 
log. Here, we have considered a 50 GeV dark matter 
particle that annihilates to bb with a cross section of 
av ^3x 10-26 cm^/s. 

For the case of a 50 GeV dark matter particle that 
annihilates to bb with av = 3 x IQ-^^ cm^/s, we esti- 
mate that FGST is expected to detect 50 or more photons 
with energy >1 GeV from only a few point-like subhalos 
in a year of observation. This result, however, depends 
strongly on the mass, annihilation cross section, and an- 
nihilation channels of the dark matter particle. If we 
increase the cross section to av = 10"^^ cm^/s (approxi- 
mately the upper limit allowed by observations of dwarf 
spheroidal galaxies [16]), for example, we expect FGST 
to detect approximately 16 subhalos, most of which with 
masses between 10^ and 10^ solar masses, prior to mass 
losses (we consider only halos up to 10^ solar masses in 
our calculations). In Fig. [2j we show the number of sub- 
halos detectable by FGST as a function of annihilation 
cross section, for selected choices of the dark matter's 
mass and dominant annihilation channel. 

At this point, we would like to comment on the un- 
certainties involved in estimating the number of subhalos 
that could potentially appear within FGST's point source 
catalog. The details of the subhalo profiles, concentra- 
tions, and the overall fraction of dark matter mass con- 
tained in subhalos can each significantly impact the num- 
ber of subhalos that are observable by FGST. In Table |lj 



we give the number of point-like subhalos observable by 
FGST for a variety of possible astrophysical assumptions, 
each for the case of a 50 GeV dark matter particle that 
annihilates to bb with av — 10^^^ cm'^/s. By varying the 
inner slope of the subhalos' density profile and the frac- 
tion of mass that is lost from the original subhalo (ML), 
we find that the number of observable subhalos can vary 
by a factor of several in either direction from the results 
shown in Fig. [2] 

If we consider the effects of substructure within subha- 
los themselves, we find that extrapolating the Bullock et 
al. concentration estimates down to a minimum mass of 
1Q~^Mq (IO-^Mq) can further enhance the number of 
observable subhalos by a factor of 2.4 (5.0). This range 
of estimates and uncertainties is consistent with those 
found, for example, in Ref. il. Furthermore, halo-to- 
halo variation in the concentration and shape of dark 
matter subhalos is expected to increase the number of 
observable subhalos by a factor of about two [T]. For 
this reason, we consider the results shown in Fig. [2] and 
Table H] to be conservative. 



III. THE FERMI POINT SOURCE CATALOG 

In this section, we turn our attention to the FGST First 
Source Catalog [2], which identifies point sources with 
greater than ~ 4cr significance within the data collected 
between August 2008 and July 2009, over the 100 MeV 
to 100 GeV energy range. Within the context of dark 
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FIG. 2: The number of point-like subhalos potentially detectable by the Fermi Gamma Ray Space Telescope as a function of 
dark matter's annihilation cross section. To qualify as detectable and point-like, we require a subhalo to produce more than 
50 events above 1 GeV per year at FGST, and require that more than 95% of those photons be concentrated within a radius 
of 2° (solid) or 1° (dotted). We show results for two dark matter masses (50 and 500 GeV), and two annihilation channels {bb 
and T^T~). 



Model/Parameters 


Detectable, Point-Like Subhalos 


Default (7 = 1.2, ML=99%) 


16.8 


7 = 1.0, ML=99% 


1.94 


7 = 1.2, ML=90% 


46.4 


7 = 1.0, ML=90% 


8.57 


Einasto a = 0.17, ML=99% 


6.50 


Einasto a = 0.17, ML=90% 


13.4 


Default w/ sub-subhalos 




M^i„ = 10-'^ Mq 


40.9 


M^in = IO-^Mq 


83.9 



TABLE I: The impact of various astrophysical assumptions 
on the number of point-like subhalos observable by FGST. 
Results are shown for a 50 GeV dark matter particle that 
annihilates with av = lO"'^^ cm^/s to 66. Note that halo-to- 
halo variation in the concentration parameter is expected to 
further increase the number of detectable point-like halos by 
a factor of ~2. See text for more details. 



matter subhalos, we are interested in sources that are 
non-variable and are not associated with known astro- 
physical objects. Furthermore, we only consider point 
sources that are more than 10° away from the Galactic 
Plane (that is, |&| > 10°), as it is expected that the re- 
gion close to the plane will contain the majority of bary- 
onic gamma ray sources (pulsars, supernova remnants. 
X-ray binaries, etc.). We also only consider sources that 
have been detected at greater than 5cr significance. Once 



these criteria have been applied, there remain 368 ob- 
jects in the FGST First Source Gatalog that qualify as 
potential dark matter subhalo candidates (from a total 
of 1451 sources). Of these, 41 have been detected with 
greater than lOcr significance. 

For each of the 368 subhalo candidate point sources, 
we compared their observed spectrum to the spectrum 
predicted from annihilating dark matter. To determine 
the observed spectrum from a given source, we took 
the gamma rays detected from a signal region within 
an angle equal to the 95% point spread function from 
the identified center of the object, and subtracted the 
average spectrum over an annulus between 3° and 5° 
around the source. Throughout our analysis, the spectra 
were binned in 34 logarithmic bins distributed between 
120 MeV and 251 GeV. For many of the point sources, 
some of the bins, especially at high energies, contained 
no events. Note that, as the point spread function of 
FGST varies considerably with energy, the angular size 
of the signal region is smaller for higher energy bins [llj . 
Events with a zenith angle greater than 105° were ex- 
cluded from our analysis, as this region is contaminated 
by gamma rays from the limb of the Earth 2J. 

Using the publicly available program PYTHIA jl2j . 
as implemented within DARKSUSY [T3], we calculated 
the spectrum of gamma rays from dark matter annihi- 
lation for 60 different masses between 10 GeV and 10 
TeV. Spectra were calculated for six annihilation chan- 
nels: bb, e+e~, fi^ ^~ , t^t~ , 'democratic leptons' (equal 
numbers of annihilations to e'^e^, iJ-^fJ.~ 
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FIG. 3: The gamma ray spectrum per dark matter annihilation for mum ~ 100 GeV, and five dominant annihilation channels; 
bb (solid black), W^W~ (dot-dashed orange), e"'"e~ (dotted red), fi^ fi~ (solid blue), and r"'"r~ (dashed yellow). 



and W~^W~ . For the gamma ray spectrum from dark 
matter annihilations to electrons, we calculated the final 
state radiation using the analytic formula described in 
Ref. |14j . In Fig. [s] we show the gamma ray spectrum 
per annihilation for the case of a 100 GeV dark matter 
particle, for several of the annihilation channels. 

Not surprisingly, we find that many of the 368 sub- 
halo candidates in the Fermi First Source Catalog can 
be well fit by a spectrum from annihilating dark matter. 
In Fig. |4l we show a few examples of such fits. Although 
the dark matter fits for each of these sources provides a 
good overall x^, one can certainly not conclude from this 
information that any of these sources are a dark matter 
subhalo. If, on the other hand, it could be shown that 
a large number of subhalo candidate sources in the cata- 
log could be well fit by a specific dark matter candidate 
(meaning a specific mass and combination of annihila- 
tion channels), then perhaps a case for a dark matter 
interpretation could be made. 

With this goal in mind, wc plot in Fig. [5] the num- 
ber of sources in the catalog that are well fit, mean- 
ing x^/D.o.F.< 1 (red) or 1.2 (black), by annihilating 
dark matter of a given mass and annihilation channel. 
Note that each point source can (and indeed does) ap- 
pear many times in each histogram. 

Examining these histograms, we notice a number of 
interesting features. Firstly, a large number of good fits 
are found for large dark matter masses. This is simply 
the result of few very high energy photons in the ob- 
served spectra of the sources. A very heavy dark matter 
particle, regardless of annihilation channel, can provide 
the required handful of ~ 100 — 200 GeV gamma rays, 
while negligibly contributing to the low energy flux. A 
class of astrophysical gamma ray sources with a very hard 
spectrum, however, could also fit the observed spectra of 
these sources. 



We also notice several bump-like features in the dis- 
tributions shown in Fig. [5] In particular, we observe 
features near 100 - 200 GeV for the bb, W+W' , and 
democratic lepton channels, and near 500 GeV for the 
T+r~ channel. We will return to possible interpretations 
of these features in Sec. El 



IV. CONSTRAINING THE DARK MATTER 
ANNIHILATION CROSS SECTION 

In this section, we use the results shown in Fig. [5] to 
determine the maximum number of dark matter subhalos 
that could plausibly be present within the Fermi First 
Source Catalog, and use that information to calculate 
an upper limit on the dark matter's annihilation cross 
section. 

To begin, we used a Monte Carlo to estimate what 
a signal from dark matter would potentially look like 
within the histograms of Fig. [5] In Fig.|6]we show, for the 
case of a 180 GeV dark matter particle that annihilates 
to W^W~ , the fraction of subhalos that are predicted to 
be fit well by an annihilating dark matter hypothesis, as 
a function of the dark matter's mass. Here we have con- 
sidered point-like subhalos that produce 50 events above 
1 GeV. As expected, the fraction with a good fit peaks 
at moM = 180 GeV, but with considerable width to the 
distribution. In the left frame, backgrounds were ne- 
glected, while in the right frame we adopted a diffuse 
gamma ray background of the form dN-^jdE^ cx E~'^ 
below 2 GeV, and dN^/dE^ oc E''^-^ above 2 GeV, nor- 
malized such that it produces 40 events above 1 GeV per 
subhalo. The inclusion of this background broadens the 
range of masses that can fit the simulated data signifi- 
cantly. At high galactic latitudes (|6| > 60°), where the 
galactic diffuse emission is lowest (~20 events per year 
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FIG. 4: Examples of annihilating dark matter fit to the spectra of objects in the Fermi First Source Catalog. The black data 
points and error bars denote the spectrum observed by FGST, while the dotted green line is the spectrum from dark matter 
annihilations, and the blue line is the sum of the background and the annihilation spectrum. See text for more details. 



per square degree above 1 GeV ) , this distribution resem- 
bles that shown in the left frame. Over most of the sky 
(20° < |6| < 60°), however, it will more closely resem- 
ble that shown in the right frame. For sources closer to 
the Galactic Plane, the distribution can be even broader. 
Note that brighter subhalos will also provide a distribu- 
tion that is more strongly peaked at the dark matter's 
mass, although we expect most of the observable subha- 
los to be not much brighter than the threshold for their 



detection. 

Using these distributions, we then proceeded to cal- 
culate how many subhalos could possibly be contained 
within the histograms shown in Fig. [Sj At a minimum, 
we could simply conclude that there there are no more 
subhalos in the sample than provide good fits. For ex- 
ample, if we consider a 400 GeV dark matter particle 
annihilating to W~^W~ , we see from Fig. [s] that about 
72 (130) subhalo candidate sources in the Fermi First 
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FIG. 5: The number of dark matter subhalo candidate point sources in the Fermi First Source Catalog that can be well fit by 
annihilating dark matter with a given mass and dominant annihilation channel. The red and black histograms denote fits with 
X^/D.o.F.< 1.0 and < 1.2, respectively. 



Source Catalog that provide a fit better than /DoF 
< 1.0 (1.2). From this, we can robustly conclude that 
fewer than approximately 145 subhalos are present in the 
catalog. 

However, if we also consider the shapes of the distri- 
butions in Fig. [5j we can in some cases more stringently 
constrain the number of subhalos contained in the cata- 



log. In particular, if a very large number of subhalos were 
actually present within one of the histograms shown in 
Fig. [5] then subtracting that population from the dis- 
tribution observed would leave a histogram with a deep, 
depression-like feature. While it is possible that the dark 
matter signal distribution (such as that shown in Fig. [6]) 
is exactly balanced by a depression in the distribution 
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FIG. 6: The fraction of simulated subhalos (for the case example of mDM annihilating to W^^Vy~, each producing 50 events 
above 1 GeV) that provide a good fit to a dark matter hypothesis, as a function of the dark matter's mass. In the left frame, 
no backgrounds are included. In the right frame, we have include a background of the form dN^/dE-y oc below 2 GeV, and 
dNry/dEj oc E~^'^ above 2 GeV, normalized such that it produces 40 events per subhalo above 1 GeV (appropriate for diffuse 
galactic emission at latitudes of 20° < |6j < 60°). 
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FIG. 7: The number of dark matter subhalo candidate point sources in the Fermi First Source Catalog that can be well fit by a 
dark matter particle annihilating to W^W~ , as a function of the dark matter particle's mass (solid). Also shown is the result 
if a population of 25 (dashed) or 50 (dotted) subhalos is subtracted from the distribution, assuming a 400 GeV dark matter 
mass. When subtracting approximately 35 or more subhalos from the distribution, a depression-like feature begins to appear 
in the distribution, indication that an oversubtraction is taking place. 



from astrophysical sources alone, we consider this unlik- 
ley and fine tuned. 

In Fig. [7j we demonstrate how we can use this ap- 
proach to more stringently constrain the number of sub- 
halos in the Fermi First Source Catalog. Here, we plot 
the histogram for the W^W" channel, subtracting from 
it a population of 0, 25, or 50 subhalos (consisting of 
dark matter particles with todm — 400 GeV). By requir- 
ing that the subtraction of a subhalo population does not 
lead to a statistically significant depression-like feature in 
the histograms, we can often place a more stringent limit 



the number of subhalos present in the Fermi First Source 
Catalog. In the case of a 400 GeV dark matter particle, 
annihilating to W~^W~ , for example, we conclude that 
less than approximately 35 subhalos are present in the 
catalog. 

We have repeated this procedure for each dark matter 
mass and annihilation channel and determined in each 
case the maximum number of dark matter subhalos pos- 
sibly contained within the Fermi First Source Catalog. 
From the maximum number of subhalos, we then cal- 
culated an upper limit on the dark matter's annihila- 
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FIG. 8: The 95% confidence level upper limits on the dark matter's annihilation cross section from studies of the Fermi First 
Source Catalog. From top-to-bottom, the curves denote dark matter annihilating to (dot-dashed red), e+e~ (dot-dashed 

black), democratic leptons (dotted red), t^t~ (dashed black), W^W~ (solid black), and bh (dotted black). We have used the 
default assumptions as discussed in Sec. |II](7 = 1.2, 99% mass loss, no sub-subhalos). See text for more details. 



tion cross section. In Fig. [8j we show the upper limit 
on the annihilation cross section as a function of mass, 
and for various annihilation channels. Here we have used 
our default assumptions as discussed in Sec.|ll (7 ~ 1.2, 
99% mass loss, no sub-subhalos), and assumed that the 
Fermi First Source Catalog is approximately complete 
for sources that produce more than ~50 events per year 
above 1 GeV. Although effects such as confusion be- 
tween multiple point sources are expected to cause some 
sources to not appear in the FGST point source catalog, 
this is thought to impact less than 10% of sources with 
|6| > 10° 0. 

Of course, if we adopt different assumptions pertaining 
to the dark matter subhalo profiles, concentrations, and 
mass losses, we can arrive at somewhat different conclu- 
sions. Erring on the optimistic side, if we use our default 
profile shape, but with only 90% mass loss, and include 
the effects of halo-to-halo variations in subhalo concen- 
trations, our predicted number of observable subhalos 
increases by a factor of ~5-6, corresponding to a limit 
on the annihilation cross section that is more stringent 
by a factor of approximately 3 (or even more stringent 
if subhalos have their own substructure). On the other 
hand, if we conservatively adopt a less cusped profile 
shape (7 = 1.0), 99% mass loss, and neglect halo-to-halo 
variations, we arrive at a limit that is about 4 times less 
stringent than that shown in Fig. [8j 

Comparing these results to other constraints that have 



been placed on the dark matter annihilation cross sec- 
tion, we find that our constraint is comparable to, or in 
some cases more stringent than, those found previously. 
In particular, the FGST collaboration has published con- 
straints from observations of 14 dwarf spheroidal galax- 
ies (Draco and Ursa Minor providing the most stringent 
constraints) [16]. For a 30 GeV dark matter particle an- 
nihilating to bb, for example, that study concludes that 
CTW < 7 X 10~^^ cm'^/s, which is comparable to our re- 
sult shown in Fig. [8] For heavier dark matter masses, 
however, the limits from dwarf spheroidals are somewhat 
more stringent (by a factor of ~ 5 at 1 TeV, for exam- 
ple). For masses below a few hundred GeV, our results 
also provide a more stringent limit than those from obser- 
vations of galactic and extragalactic diffuse emission [17] . 



V. POSSIBLE HINTS OF DARK MATTER IN 
THE FERMI FIRST SOURCE CATALOG? 

In this section, we turn our attention to some of the 
features observed in the histograms of Fig. [5] and ex- 
plore the possibility that perhaps ^10-50 of the sources 
in the Fermi First Source Catalog might be dark matter 
subhalos. 

In particular, we notice in Fig.[5]potentially interesting 
bump-like features at todm ~ 500 GeV, for annihilations 
to T+T", and at ttidm ~ 150 — 200 GeV, for annihila- 
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FIG. 9: The number of dark matter subhalo candidate point sources in the Fermi First Source Catalog that can be well fit by 
a dark matter particle annihilating to t'^t~ (left) or W^W~ (right), as a function of the dark matter particle's mass (solid). 
Also shown in each frame is the result if a population of 25 (dashed) or 50 (dotted) subhalos is subtracted from the distribution, 
assuming a 500 GeV or 180 GeV dark matter mass in the left and right frames, respectively. When subtracting approximately 
30 subhalos from either distribution, the observed bump-like features are flattened. The observed bump-like feature could thus 
be explained if approximately 30 of the sources within the Fermi First Source Catalog are dark matter subhalos. 



tions to W^W^ , bb, or democratic leptons. Could these 
features be the result of a population of dark matter sub- 
halos within the FGST point source catalog? To assess 
this question, we subtracted from the observed distribu- 
tion of sources that which would be predicted (by our 
Monte Carlo) from a dark matter subhalo population. 
In Fig. [9j we show example results of this subtraction. 
For the case of a 500 GeV dark matter particle annihi- 
lating to r+T~, we find that the observed bump can be 
removed, flattening the overall distribution, if a popula- 
tion of approximately 30 subhalos is present within the 
catalog. A similar number of subhalos with todm = 180 
GeV, annihilating to either W^W^ or bb, would also 
also largely removed the features seen in their respective 
channels. 

In the mDM = 500 GeV, r+r" case, in order to pro- 
duce the required number of observable subhalos, the 
dark matter must have a very large annihilation cross 
section (crw ^ 6 x 10^^'^ cm'^/s). Intriguingly, in order for 
such a dark matter particle to provide the excess positron 
fraction reported by PAMELA 18J, it must have a cross 
section of approximately cru ~ 2 x 10"^'^ [12]. The factor 
of ^3 disparity between these two values could easily be 
accounted for if, for example, subhalos contain significant 
substructure, or less mass loss has taken place than we 
had assumed (see Table |l]). 

In the case of todm = 180 GeV, and annihilations 
to W^W'^ , we find again that a population of approxi- 
mately 30 subhalos within the Fermi First Source Catalog 
can flatten the observed distribution. This corresponds 
to an annihilation cross section of ^ (1—2) x 10^^^ cm'^/s. 
This is comparable to the value predicted for a wino-like 
neutralino, such as is found in models of anomaly medi- 
ated supersymmetry breaking [20] . A similar cross sec- 
tion would also be found for the case of annihilations to 



bb. We do not consider the democratic lepton case further 
as this appears to conflict with the cosmic ray electron 
spectrum as measured by FGST |5T]. 

We emphasize that none these features can, at this 
time, be used to confidently support a claim that a popu- 
lation of dark matter subhalos is present within the Fermi 
First Source Catalog. It is the case, however, that the 
presence of such a population could explain some of the 
features observed in the catalog. If either of these ob- 
served features does, in fact, result from the presence of 
a subhalo population, the feature is be expected to be- 
come more prominent as the FGST accumulates more 
data, and the subhalos' spectra become more precisely 
measured. 



VI. SUMMARY AND CONCLUSIONS 

In this article, we have studied the recently published 
Fermi Gamma Ray Space Telescope's First Source Cata- 
log, and considered whether a significant number of the 
sources in this catalog might be dark matter subhalos. 
For a typical thermal dark matter candidate with a mass 
of 50 GeV, and reasonable astrophysical assumptions, 
one predicts that a few relatively large 10"^ — lO^M©) 
and nearby (~ 0.01 — 10 kpc) subhalos would be bright 
enough to appear within such a catalog. It is unlikely, 
however, that such a small number of sources could 
be identified from among the hundreds of unidentified 
sources in the Fermi First Source Catalog. If the dark 
matter's annihilation cross section were larger than that 
typically predicted for a thermal relic, however, it is pos- 
sible that a larger number of subhalos could be detectable 
by gamma ray telescopes. 

From among Fermi's First Source Catalog, we studied 
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368 sources as dark matter subhalo candidates, each of 
which has been detected with high significance (> 5cr), is 
located away from the Galactic Plane {\b\ > 10°), and is 
not associated with any source observed at other wave- 
lengths. Although we find that the spectrum of many 
of these sources could be well fit by that predicted from 
annihilating dark matter, this is most likely the result of 
the relatively large error bars involved, and on the wide 
range of dark matter masses and annihilation channels 
we have considered. 

By studying the distribution of source spectra in the 
First Source Catalog, we derived an upper limit on the 
annihilation cross section of dark matter, as a function 
of its mass and dominant annihilation channel. For dark 
matter particles with relatively low masses (less than a 
few hundred GeV) , we find constraints that are compara- 
bly stringent to those derived from observations of dwarf 
spheroidal galaxies. For example, in the case of a 50 GeV 
dark matter particle that annihilates to bb, we conclude 



that av <7 X lO'^^ cm^/s. 

Lastly, we note the appearance of bump-like features 
in the distribution of spectra observed in the Fermi First 
Source Catalog. Such features could be explained (ie. 
flattened) if the dark matter takes the form of a 500 GeV 
particle with an annihilation cross section to r+r" of 
av ^ 6x 10"^'^ cm'^/s, or a 180 GeV particle with an an- 
nihilation cross section to W~^W~ of ctw ~ (1 — 2) x 10"^'* 
cm^/s (such as is predicted for a wino-like neutralino). 
These scenarios could also potentially account for the 
cosmic positron excess observed by PAMELA. 
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